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Abstract 

We compute the 1/N correction to the location of the previously found first-order 
phase transition in the Gross-Neveu model at a chemical potential fj, = /i c = -^j? 7 ^ where 
m is the fermion mass. We employ an expression for the free energy /(/i) given by the 
thermodynamic Bethe ansatz under the approximation that the fundamental fermions 
dominate the ground state, and combine it with the effective potential evaluated at zero 
chemical potential. Our result is \x c = 77= [1 — ^^]- 
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1. Introduction 

In a previous paper |TJ, we studied the Gross- Neveu model 0] in the presence of a 
chemical potential \i. We computed the effective potential analytically to leading order in 
1/N, and we found a first-order phase transition at 

1 

H = = -j=m. (1) 

where m is the fermion mass. For fx > fx c the discrete chiral symmetry is restored and 
the fermion becomes massless. The earlier references on the phase transition of the Gross- 
Neveu model may be found in ref. Q. 

It is known Q || that the N = 2 Gross-Neveu model provides an approximate 
description of trans-polyacetylene |3|, and we used our results to explain an observed 
phase transition in polyacetylene as a function of doping concentration (which is related 
to the chemical potential) 0]. We found agreement with the order, the nature, and the 
location of the phase transition, and we concluded that the Gross-Neveu model did indeed 
capture some of the essential physics of trans-polyacetylene. 

One aspect of this work that requires further amplification, however, is the fact that 
polyacetylene corresponds to the N = 2 Gross-Neveu model, whereas our results were 
derived only to leading order in 1/N. This raises the possibility that the next leading term 
in the 1/N expansion might provide a quantitatively significant correction, which would 
then render the excellent agreement between theory and experiment found in leading order 
merely fortuitous. 

In this paper we address the question of the order 1/N correction to our previous 
computation. The straightforward way to do this would be to evaluate the Feynman 
graphs that contribute to order 1/N, incorporating, as we did in leading order, the effect 
of \x into a modification of the fermion propagator. However, the evaluation of the 1/N 
correction to the effective potential even when \i = is a rather involved affair, and the 
extension to non-zero \i appears to be prohibitively complicated. Instead, we shall evolve 
a technique based on the thermodynamic Bethe ansatz (TBA) || 0. 

In the next section, we shall explain this technique and illustrate it by rederiving the 
results of our previous paper. Then in sections 3 and 4, we shall extend the computation 
to next order in 1/N; for simplicity, we shall find it convenient to make a further approxi- 
mation in which we restrict the Gross-Neveu S- matrix to a single component. This will be 
explained more fully below. In section 5 we shall summarize our results and draw some 
conclusions. 
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2. The TBA in Leading Order 

The TBA has been much used recently to discuss the properties of exactly integrable 
models H, jlCR. In this paper, somewhat unconventionally, we shall actually use it to 



determine the thermodynamic properties of a particular system. The TBA should be very 
useful technique in discussing thermodynamics of the Gross-Neveu model because all the 
S- matrix elements are known [JTJ |T^] ||13| . 

The basic quantity of interest is the " dressed" excitation energy e(6>) that is required 
to promote a single particle from the Dirac sea to an energy mcoshO; e{9) includes the 
effects of the rearrangement of the vacuum after the particle has been excited. At finite 
temperature, the equation for e(9) is non-linear, but at zero temperature (which is the case 
of interest to us) it reduces to the linear equation [TJ 



e(9) = li - mcosh9 + [ d9' K{9 - 9')e{9'). (2) 

J-B 

Here the kernel K is the logarithmic derivative of the S-matrix: 
and B is determined by the condition 

e(±B) = 0. (4) 

Effectively this condition determines B as a function of the chemical potential \x. 

Once e is known, the free energy / of the system can be determined as a function of 
\i up to an additive constant: 

/B 
d9e(9)cosh9. (5) 
-B 



As shown in ref. |Tl[ , the expansion of K in powers of 1/N begins in order 1/N, so to leading 



order one ignores the integral on the right-hand side of eq. (2), and one has immediately 

e(9) = fi — mcosh9 (6) 



and hence 



coshB = fx/m. (7) 
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This is possible only if jj, > m; if jj, < m, no solution for e exists. 
^From eq. (5) we find 



Til 1 

/(aO - /(0) =^ 2 - m 2 )( — )(B - -«nfe2B) 

27r 2 2 (8) 

=9(/j 2 — m 2 )(——)[ln(fi/m + a/ [i 2 jm 2 — 1) — fi/m-s/ fj^/rri 2 — 1]. 
On dimensional grounds, /(0) = —bm 2 where 6 is a dimensionless constant, i.e. 



/(//) = —bm 2 + 9(fx 2 — ?™ 2 ) — — [ln(fi/m + a/ \i 2 jm 2 — 1) — \i j m\/ \i 2 j m 2 — 1]. (9) 

At first sight, there appears to be no way that this expression can be used to provide 
evidence for a phase transition at fi = m/y/2, since /(//) is absolutely flat for \i < m. We 
observe, however, that if a massless phase exists, then its free energy fo(fi) is given by eq. 
(9) with m=0: 

foQj.) = -v 2 /2n. (10) 
Provided that < b < 1/2%, the curves /(//) and fo([i) will intersect at 

li 2 /m 2 = 2tt6, (11) 

with fo(fi) > f(n) for (///m) < 27r6. One therefore predicts a first-order phase transition 
from the massive phase to the massless one at \xjm = V2nb. 

The one piece of information the TBA does not give us is the value of b. For this, 
we need the effective potential computation, but only at \x = 0. Thus the TBA plus the 
effective potential at \x = allows us to determine the phase structure of the theory for 
H > 0. 

To make contact with previous work, we note that the free energy as given above 
corresponds to the value of the effective potential V at its minimum a m i n , and the mass 
m is equal to gamin where g is the GN coupling constant. We can verify this explicitly by 
appealing to eqn. (20) of ref. Q, from which we obtain 
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V e f f (a) = ^ + ^[6(a 2 - 7 2 )(/n(a 2 /a 2 ) - 3) + 0( 7 2 - a 2 )(2/n 7 + ^ - 3)] 



^Vt^^t 2 -^ 2 )- 

(12) 

This expression comes equipped with the following comments: 

(i) We have divided V e f / of ref. [1] by AT. This ensures a finite limit as N — ► oo (with 
g 2 N and m fixed) and also is the correct normalization to correspond to /(//), as we shall 
see explicitly below; 

(ii) We have omitted the last term of eq. (20) of ref. [1] , since it is precisely cancelled 
by another contribution that was left out of eq. (20), called V in eq. (14) of ref. 0. Thus 
V e ff above includes all the relevant contributions; 

(iii) The parameter 7 above is fiy/N/n. 

(iv) We have set the coupling A = g 2 N that appears in eq. (20) of ref. |l|equal to tv. 
This is permissible because physics depends only on 

CTrmn = CT^ 1 ^'^ (13) 

where ao is an arbitrary renormalization point. The choice A = n is convenient because 
then a mi n = ctq. We see that 



Wff0) = 4 + w _ .ai^nfS, _ 2vp, . (14) 

2 2 

Putting m 2 = g 2 cr 2 = ^f-, and jj? = we have 



m 2 2 2\ m2 /U + a//U 2 — m 2 \i In 2 



Veff = ~ + - m')^-[ln^-^ — -rL,r_-l] (15) 

47T 2n m m V m 1 

which tallies exactly with the expression for /(//), as given by the TBA, provided b = l/4n 
(and hence \i c = m/y/2). 

Of course, as pointed out above, we did not need the full V e ff, including the effects 
of n, to determine b. It would have sufficed to know the original V e ff(a), computed by 
Gross and Neveu 0: 

Vgn(*) = ^ + ^[ln(a 2 /a 2 ) - 3] (16) 
from which Vqn{g = °"o = \ v m ) = ~~ Tn f°ll° ws immediately as well. 
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3. Next to Leading Order 

In next order, we of course do not have the explicit form of V e f / as a function of a and 



fi. For \x = 0, the 1/N correction was computed long ago by Schonfeld |15j and by Root 
||l(j|| . Their results will allow us to extract the 1/N correction to the value b = found 
in the previous section. The other ingredient we shall need is the form of f(ji) obtained 
from the TBA, especially the massless limit /o (/•*)• We shall consider each of these in turn. 

For our purposes, it is convenient to follow the formulation of the 1/N correction AV 
to the Gross-Neveu model given by Schonfeld. By summing Feynman diagrams, he obtains 
an unrenormalized expression. 

AV un = [ dx\ln[-A + ln(a/a )}-lnln(-^)\ (17) 
Jo I 2 9 ) 

where we have set A = tv in Schonfeld's expression, and where A = Z ln(%^j), and 
Z = [1 + 4g J 7 j 1 / 2 . AV un is divergent, and it also has an imaginary part. As Schonfeld 
explains, the correct physical prescription, at least in the range of a we are interested 
in (cr ~ ctq) is to take the real part. The divergences are removed by a renormalization 
prescription which amounts to the conditions 

d 2 

ReAV \ a=!J = —(ReAV) \ a=ao = . (18) 
Now in leading order, a m i n = cxo (since A = tt), and so we must have 

<T min = (ro + jj(ri, with V = V GN + ReAV (19) 

and 

V(a min ) = V GN (a + -^£7i ) + ReAV(a + -^<ri) . (20) 
But V GN (ao) = 0, and AV is already order 1/N, and so to the order we are working, 

V(a min ) =V GN {a ) + ReAV {(j ) 

=Vgn(&o) because of eq. (18) 



(21) 



4N 

4N ^ N^ 1 ^j^\®min j^&min&l, 



The quantity o\ is determined from the condition 



V'(a + ^a 1 ) = (22) 

and is given by 

£ = -ReAV'(ao) (23) 

where use has been made of the renormalization condition V" (erg) = 1. 

It can easily be seen from Schonfeld's paper that in terms of the unrenormalized 
expression AV un defined above, 

ReAV'(ao) = ReAV ,un {a Q ) - ^-ReAV un {a Q ) - \a Q AV" un {a Q ) (24) 

oo"o 6 

where the last two terms on the rhs are the counterterms necessary to implement the 
conditions (18) and thereby render the expression finite. 

The remaining task is the purely numerical one of using the integral expression (17) 
to evaluate the right-hand side of eq. (24) . One finds 

°"i ,1.06. . . 

N = < W )ff ° (26) 

and hence 

. . 1 2 1 ,1.06, 2 /nr\ 

^(%m) = -^^ m ,n+ ^2 (y l^mrn ( 26 ) 

where in the second term we have used a m i n = <jq + $(l/iV). 
Hence V min = -bm 2 where b = ^[1 - |^]. 
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4. The TBA in Next to Leading Order 

Let us formulate the TBA in the framework of 1/N expansion. While we actually deal 
only with the next to leading order in this paper, we present the formulas in an organized 
way so that the structure of the higher orders in 1/N becomes transparent. We expand 
various quantities introduced in section 2 in powers of 1/N: 

g(0) = g o (0) + l e ~ 1 (0) + ... 

M-f{Q) = M + ^gM + --- (27) 

One might think that the boundary of the integral B in (2) determined by (4) receives a 
1/N correction and that it is therefore incumbent on us to expand B in powers of 1/N. 
However, by examining a model in which the TBA equation can be solved exactly, we learn 
that it is more appropriate to keep B fixed and change \i as iV is varied. We will treat this 
in the Appendix. Thus we expand the chemical potential 

f ,(B)= f i (B) + ^ f i 1 (B) + ... . (28) 
Ho is determined by the leading order equation 

e (B)=0 . (29) 

For clarity we recapitulate the leading order result obtained in section 2 in the present 
notation: 

= Ho ~~ mcosh9 

m 2 1 ( 30 ) 
go{H) =0(v 2 o-m 2 ) — (B--sinh2B) . 

Ztv Z 



Following Forgacs, Niedermayer and Weisz [14], we make the ansatz that the ground state 



of the Gross- Neveu model at finite density and zero temperature is dominated at large- N 
by the fundamental fermions. This is a reasonable assumption because the solitons and 
their bound states are more massive than the fundamental fermions at large N. This 
ansatz allows us to ignore all the S-matrix elements but the symmetric combination of the 



fundamental fermion scattering amplitudes, as shown in [14|. The approximation is in fact 
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implicit in (2); without the ansatz K becomes a matrix which acts on a vector e{9) in the 
particle species basis. 



The consistency of the ansatz is checked by Forgacs, et al., in refs. |14[] and Wn - The 
perfect matching of the large- N and the TBA results in the leading order, as discussed in 
the preceding section, lends further support to this viewpoint. 

With this ansatz, the kernel K in the TBA equation (3) in leading non-vanishing order 
has the form 

K M = -S<? " 55S>- (31) 

Using (2) it is straightforward to compute the 1/N correction of the dressed particle energy 
e(6>) and the free energy f(p) — /(0). They read 



ei =fi! + msinh 9[Chi{9 + B) - Chi(6 - B) - In 



sinh{9 + B) 



sinh(6 - B) 

- mcosh 6[Shi(9 + B)- Shi{9 - B) - 2B] (32) 



gi (jj) = - 0(p 2 - m 2 ) — [B 2 + sinh 2 B - BShi(2B)} 



TV 



where Shi(x) and Chi(x) are the hyperbolic integral functions defined by 



, rX sink t 
bni[x) = at 



o 



Chi(x) = 7 + Inx + dt 



t 

cosh t — 1 



(33) 



Here 7 is Euler's constant. The result for <7i(//) in (32) reproduces that of Forgacs, et al. 



The boundary condition e(B) = with €q(B) = 0, which is already met in the leading 
order, determines fio and //].: 

uo = mcoshB 

(34) 

fH = mcoshB [Shi(2B) - 2B ] - msinhB [Chi(2B) - In {sinh2B} - 7] 

To implement the strategy explained in section 2 we need the expression for the free energy 
in the massless limit, which corresponds to B —>■ 00. Using the asymptotic form of the 
hyperbolic integral functions 
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e x 



Shi(x) — > x ^oo ^ 



(35) 



one can easily show that 



Chi(x) — >£— >oo — 



lim . gx{n) = 

m— »0 

lim Hi = //o(7 — ^2) 

m— »0 



(36) 



The latter implies that 



^ = M1 + ^^) • 07) 



Therefore, we obtain 



to[/(,)-/(0)] = -| 

2tt V iV ; 

We deduce by combining (26) and (38), that the phase transition takes place at 

£ 0232 ^2 

27T iV 47T 3iV V ; 



or 



If we estimate (40) at N = 2 we have 

/*= ^(1-0.23) (41) 

which implies about a 20% correction to the critical chemical potential due to the next to 
leading order in 1/N. 
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5. Conclusions 

In our previous paper, in leading order in the 1/N expansion, we found a first-order 
phase transition in the Gross-Neveu model as a function of chemical potential, occurring 
at ix c = T^m. 

In this paper, making use of the 1/N corrections to the fx = effective potential, and 
also of the thermodynamic Bethe ansatz, we have shown that the above result receives a 
correction of the form 

"-7f mll -^l ' (42) 

In view of the result, it is reassuring to note that the correction we have found acts to 
lower the value of fi c , i.e. to make it easier for the system to undergo the phase transition. 
The validity of our picture might begin to run into trouble if fi c were to become larger 
and were to approach m in magnitude, because for \i > m /(/x) decreases rapidly from its 
constant value of /(0), and the question of whether f{pt) and /o(a0 intersect at all would 
then become more problematical. 

When we evaluate our correction for N = 2 we find a 20% effect. Since the value 
of the doping concentration y is, in this model, directly proportional to \i, this implies 
that y changes from .064 to perhaps .050, which is in somewhat less good agreement with 
published results, but still well within the experimental limits. 

In obtaining the above result, we have made one further approximation: following 
Forgacs, et al., we have restricted the S-matrix appearing in the TBA equation to a single 
component, the scattering between fundamental fermions. We believe that it is a good 
approximation at large-N because all other excitations but fundamental fermions have 
masses of the order of N. In fact, the result by Forgacs et al. obtained under the same 
approximation indicates a perfect matching between the TBA and the perturbative cal- 
culations even at large values of chemical potential. Nonetheless, we should mention that 
the validity of the approximation has to be checked by more elaborate computation which 
includes all the solitons and bound states. The obstacle at present to such computation is 
a considerable amount of complexity and labor. 

Also as we noted above, we have chosen what may seem a curious strategy in perform- 
ing the 1/N expansion for the TBA. Rather than keeping the physical parameter ji fixed, 
and developing the auxiliary variable B in a series in 1/N, we have inverted this procedure. 
Of course, the two must ultimately be equivalent, but we have found, in studying a simple 
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exactly solvable example of the TBA equation, that the latter is more transparent and 
more manageable. We amplify this point in a brief appendix. 

Another aspect we have discussed in this paper is how the density phase transition 
between massive and massless phases can be addressed within the framework of the TBA. 
The TBA itself only allows us to evaluate the difference between free energies in finite and 
zero chemical potential theories. We have proposed a new machinery of combining the 
massless limit of the TBA result with that of the computation of the effective potential 
at zero chemical potential. We have verified that it works in the leading order in the 1/N 
expansion of the Gross-Neveu model. It remains to be seen if it has wider applicability. It 
would also be a challenging task to see whether one can develop a formalism that could 
reveal this type of phase transition solely within the framework of the TBA. 
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Appendix A. 

We give a brief description of an exactly solvable model of the TBA. The model is 
defined by the 2-body S-matrix 

S ab = exp(^£e^p%p%), (Al) 

TXX 

where A is a real constant. The S-matrix satisfies unitarity and "crossing" but does not have 
polynomial boundedness, and hence may not be realized in local quantum field theories. 
We shall utilize the model as a theoretical laboratory to examine the large-N expansion 
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of the TBA. We nevertheless keep in mind the possibility that it may allow a physical 
interpretation. Properties of models without polynomial boundedness have recently been 
discussed by Khuri [JTBfl . With the definition = m(cosh9i, sinh9i), i = a,b, and the 
introduction of the relative rapidity 9 = 9 a — Of,, the kernel K takes the simple form 

K{0) = XcoshO. (A2) 
The TBA equation (2) can be readily solved as 

e(0) = [A — mcoshO (A3) 

where 

it) 

fa = = (A4) 

1 + X[\sinh2B - B] y ' 

Notice that B is determined by (4), \i — rhcoshB = 0, and we can obtain B nonperturba- 
tively in this model. The free energy of the model can be easily computed as 

1 

m - /(0) = -e{fi - m)( — ) — . (A5) 

27r A + l^sinhzB — B\ 1 

We treat the model perturbatively to gain insight for the large N expansion in the 
TBA. We expand various quantities in powers of A: 

fjL = Ho + A/ii + A 2 /U 2 + • • •• 
e(0) = e {9) + Xh{e) + \ 2 e 2 (6) +■■■. (A6) 
f(fi) - /(0) = 9o (h) + XgM + X 2 92(h) + ■■■■ 

But we keep B fixed as in our treatment of the Gross-Neveu model in section 4. Then one 
can solve the TBA equation order by order in perturbation theory to obtain the results: 

[in = fio[B - ]-sinh2B] n , 

2 J (AT) 
~e n (e) = e (8)[B--sinh2B} n , 

where 

u = mcoshB, 

(AS) 

€q(6) = Ho — mcoshO = m(coshB — coshO). 
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It can be readily seen that they reproduce the exact solution (A3) and (A4). 

On the other hand, the perturbative treatment of the model when B is expanded as a 
power series in A is far from transparent and does not appear to work, as may be discovered 
by making the attempt. 

There is a whole family of models (of which the above is the simplest member) for 
which the TBA can be solved in similar fashion. 



13 



References 



[1] A. Chodos and H. Minakata, Phys. Lett. A191 (1994) 39. 
[2] D.J. Gross and A. Neveu, Phys. Rev. D10 (1974) 3235. 

[3] R. F. Dashen, S.-K. Ma, and R. Rajaraman, Phys. Rev. Dll (1975) 1499; U. Wolff, 

Phys. Lett. B157 (1985) 303. 
[4] H. Takayama, Y. R. Lin-liu, and K. Maki, Phys. Rev. B21 (1980) 2388. 
[5] D. K. Campbell and A. R. Bishop, Nucl. Phys. B200 (1982) 297. 
[6] W. P. Su, J. R. Schrieffer, and A. J. Heeger, Phys. Rev. B22 (1980) 2099; A. J. Heeger, 

S. Kivelson, J. R. Schrieffer, and W. P. Su, Rev. Mod. Phys.60 (1988) 781. 
[7] J. Chen, T.-C. Chung, F. M. Moraes, and A. J. Heeger, Solid State Commun. 53 

(1985) 757; F. M. Moraes, J. Chen, T.-C. Chung, and A. J. Heeger, Synth. Met. 11 

(1985) 271. 

[8] C. N. Yang and C. P. Yang, J. Math. Phys. 10 (1969) 1115. 

[9] Al. B. Zamolodchikov, Nucl. Phys. B342 (1990) 695. 
[10] T. R. Klassen and E. Melzer, Nucl. Phys. B338 (1990) 485; B350 (1991) 635. 
[11] A. B. Zamolodchikov and Al. B. Zamolodchikov, Ann. Phys. (NY) 120 (1979) 253. 
[12] R Shankar and E. Witten, Nucl. Phys. B141 (1978) 349. 
[13] M. Karowski and H. J. Thun, Nucl. Phys. B190 [FS3] (1981) 61. 
[14] P. Forgacs, F. Niedermayer, and P. Weisz, Nucl. Phys. 367 (1991) 123. 
[15] J. Schonfeld, Nucl. Phys. B95 (1975) 148. 
[16] R. G. Root, Phys. Rev. Dll (1975) 831. 

[17] P. Forgacs, F. Niedermayer, and P. Weisz, Nucl. Phys. 367 (1991) 144. 

[18] N.N. Khuri, |hep-ph/9405406| , |hep-ph/9512386| (to be published in proceedings of the 



"Vlth Blois Workshop, Frontiers in Strong Interactions," June 1995). 



14 



